Molecular Conductance: Chemical Trends of Anchoring Groups 
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Combining density functional theory calculations for molecular electronic structure with a Green 
function method for electron transport, we calculate from first principles the molecular conductance 
of benzene connected to two Au leads through different anchoring atoms - S, Se, and Te. The 
relaxed atomic structure of the contact, different lead orientations, and different adsorption sites 
are fully considered. We find that the molecule-lead coupling, electron transfer, and conductance 
all depend strongly on the adsorption site, lead orientation, and local contact atomic configuration. 
For fiat contacts the conductance decreases as the atomic number of the anchoring atom increases, 
regardless of the adsorption site, lead orientation, or bias. For small bias this chemical trend is, 
however, dependent on the contact atomic configuration: an additional Au atom at the contact with 
the (111) lead changes the best anchoring atom from S to Se although for large bias the original 
chemical trend is recovered. 

PACS numbers: 73.40.Cg, 72.10.-d, 85.65.-|-h 



I. INTRODUCTION 

A critical issue in molecular electronicsiiSiii is to find 
anchoring groups and construct contact structures which 
provide both stability and high contact transparency. 
In many recent experimentSf^iSiSiS Au electrodes were 
used as leads for electronic current because of their high 
conductivity, stability, and well-defined fabrication tech- 
nique. A common way to construct a lead-molecule-lead 
(LML) system is by using a break junction, formed ei- 
ther mechanicallji^iSiL2iia or electrically i^iiiJ^ In these 
break-junction experiments, the atomic structure of the 
molecule-lead contact is unknown. Therefore, neither 
the influence of detailed atomic structure on transport 
through the molecules nor a path to improved perfor- 
mance is clear. With regard to anchoring groups, a 
S atom is commonly usec&^iSiLSiS to connect various 
benzene-like organic molecules to Au leads because of 
its good binding with Au and good stability. However, a 
recent experiment for a,w-bisacetylthio-terthiophene and 
a,w-bisacetylseleno-terthiophen9i4 showed that a Se an- 
choring atom is better than S for zero bias conductance. 
Currently, it is an open problem whether there are other 
better choices and what the system dependence is. 

In terms of theoretical studies, there are mainly 
two ab initio approaches for electron transport through 
molecules. One was developed by Lang, et ali^: the 
Kohn-Sham equation of the system is mapped into 
a Lippmann-Schwinger scattering equation which is 
solved for the scattering states self-consistently. In the 
implementationfiS*i& the jellium model was adopted for 
the two metallic electrodes of an LML system. The other 
approach fiLi24i2i22i2iiSS is based on a density functional 
theory (DFT) calculation for the molecular electronic 
structure combined with a non-equilibrium Green func- 
tion (NEGF) method for electron transport. Very close 
to the latter approach there is also an approach based 
on a self-consistent tight-binding method combined with 



the NEGF method,^i24 

Within the DFT+NEGF method itself, there are gen- 
erally two approaches describing the boundary between 
the leads and the molecule of a LML system. In one cat- 
egory, a duster geometry is adopted, either for all the 
subsystems of a LML system or for only the extended 
molecule while the leads are handled by an extended 
system method (for example, Refs.)iLi2ii2*2SiSL21 The 
other type of DFT+NEGF approach22i2ii2^ uses periodic 
boundary conditions (PBC) with large parts of the leads 
included in the extended molecule, so that the interac- 
tion between the molecule and its images will be screened 
off by the metallic lead in between. In this case all the 
subsystems can be treated on the same footing. 

Concerning the important issue of the best anchoring 
group, surprisingly few studies have been done. First, 
there is a theoretical calculationi^*i& for the benzene 
molecule using jellium leads which demonstrated that the 
best anchoring atom for contact transparency is not an S 
atom but rather a Te atom - it increases the equilibrium 
conductance by about 25 times. A potential problem 
with this calculation, however, is the neglect of the de- 
tailed electronic structure of Au as well as the contact 
atomic structure. Although in the DFT-I-NEGF method 
both molecule-lead interaction and contact structure can 
be taken into account in principle, in practice this has 
not been done: contact atomic structures were prede- 
termined in all previous calculations. Thus, to the best 
of our knowledge, the effect of relaxing atomic positions 
in the contacts has not been investigated theoretically. 
Very recently, the effect of three different molecular end 
groups has been evaluated with DFT-t-NEGF)2& though 
not the ones considered here nor with contact structure 
relaxation. 

In this paper we investigate systematically the ef- 
fects of different realistic contact atomic structures and 
re-examine the important issue of the best anchoring 
group. We use our previously developed self-consistent 
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implementation of the DFT+NEGF methodiSS In our 
method, as in the other two implementations j^SiSi PBC 
are adopted and large parts of the two metallic leads 
are included in the device region (about 40 Au atoms 
per lead) so that the molecule-lead interactions (includ- 
ing electron transfer and atomic relaxation) are fully in- 
cluded, and the electronic structure of the molecule and 
leads are treated on exactly the same footing. 

We calculate the molecular conductance of benzene 
connected to two Au leads with real atomic structure and 
finite width in order to simulate physical leads in break- 
junction experiments. The connection is made through 
three different anchoring atoms: S, Se, and Te. We con- 
sider fully the atomic relaxation of these LML systems, 
and consider different lead orientations and different ad- 
sorption sites of the anchoring groups. The fully opti- 
mized atomic structures and high level basis set for all 
the atomic species make our calculations well converged. 
It is found that the molecule-lead coupling, the electron 
transfer, and therefore the conductance depend strongly 
on the adsorption site, lead orientation, and local con- 
tact atomic configuration. It turns out that for ideally 
flat contacts the equilibrium conductance decreases with 
increasing atomic number of the anchoring group, regard- 
less of the contact atomic structure and bias. However, 
for contacts with a fluctuation of one Au atom, Se be- 
comes slightly better than S for small bias, although for 
large bias the chemical trend goes back and S becomes 
better. This shows the critical role of the real contact 
atomic structure in electron transport through molecules. 



II. COMPUTATIONAL METHODS 

The systems studied contain a benzene molecule con- 
nected to two Au leads through S, Se, or Te anchoring 
atoms. While our main focus is on the chemical trends 
among these anchoring atoms, the structure of the con- 
tact is an important issue as well. Two different Au lead 
orientations, (001) and (111), are considered. The in- 
plane size of the Au lead is set to be 2^/2 x 2\/2 for the 
(001) lead and 2x2 for (111). For checking the depen- 
dence on the in-plane lead size, we also use larger (001) 
leads, 3-\/2x3a/2, and 4-\/2x4-\/2. Furthermore, we also 
use a (001)-(4-\/2 X 4-\/2) periodic surface to model the 
(001) lead. Because of the large separation between the 
molecule and its images (larger than 12A), this model 
will be a good approximation to the lead of an infinitely 
large surface. For the (111) lead, we consider different 
adsorption sites for the anchoring atom: hollow site (h), 
bridge site (b), and top site (t). From now on we use, 
for instance, the structural label S_(001)_h to denote the 
system with S anchoring atom adsorbed at the hollow 
site of the (001) lead surface. 

The purpose here is to simulate possible experimental 
situations in break-junction experiments, in which dif- 
ferent atomic structures may occur and, indeed, global 
structural equilibrium may not be reached. In our calcu- 
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FIG. 1: Optimized atomic structures for systems with two 
S anchoring atoms. Optimized structures with two Se or Te 
anchoring atoms are similar and therefore not shown, (a) The 
Au leads are in the (001) direction, and the anchoring atom 
is located at the hollow site of the Au(OOl) surface [as shown 
in (e)]. (b)-(d) The Au leads are in the (111) direction, and 
the anchoring atom is located at the hollow, bridge, and top 
sites, respectively [as shown in (f)]. The dashed line denotes 
the interface between the device region (C) and the left or 
right lead {L or R). 



lation, the contact atomic relaxation is fully included by 
relaxing the molecule and the first two atomic layers of 
the lead surfaces, as well as the molecule-lead separation, 
while leaving the in-plane position of the anchoring atom 
fixed at the h, b, or t site. The optimized atomic struc- 
ture of the systems with an S anchoring atom are shown 
in Fig. 1 ; those for the other anchoring atoms are similar 
and therefore not shown. These fully optimized struc- 
tures are called "relaxed" structures hereafter (default 
if not specified), in contrast to "unrelaxed" structures in 
which atoms in the leads are fixed at their bulk positions. 
In the unrelaxed case, the molecule is fixed at its opti- 
mal isolated structure, the dangling bond on the S atom 
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TABLE I: Binding energies between the lead and the molecule (in eV), in terms of the total energies {Etot's) of the whole and the 
subunits of a system. For a relaxed or unrelaxed system it is determined by _Etot(load+molecule) - iJtot (lead) - Btot (molecule). 













Se- anchored 




Te- anchored 


orient. 


site 


relaxed 


unrelaxed 


relaxed unrelaxed 


relaxed unrelaxed 


(001) 


h 


3.40 


3.41 


3.28 


3.32 


3.07 


3.07 




b 


2.63 


2.50 


2.51 


2.42 


2.45 


2.38 




t 


1.60 


1.60 


1.49 


1.55 


1.66 


1.65 


(111) 


h 


4.23 


4.27 


4.10 


4.11 


3.89 


3.88 




b 


4.44 


3.86 


4.33 


3.72 


4.14 


3.55 




t 


3.55 


2.60 


3.52 


2.55 


3.38 


2.67 



is saturated by an Au atom, and the distance between 
the S atom and the Au surface is optimized. 

Because we use the bulk Au structure for the leads, 
the atomic relaxation consists of two parts: one is the 
relaxation of the bare Au lead with respect to its bulk 
structure, and the other is the relaxation of both the 
leads and the molecule induced by the molecule-lead in- 
teraction. The latter is our primary interest; in fact, our 
calculations show that the former is very small, which is 
consistent with the very small surface relaxation of un- 
reconstructed infinite Au surfaces. 

For the electronic structure calculation, we use Siesta, 
an efficient full DFT package.-^ A high level double 
zeta plus polarization basis set (DZP) is adopted for 
all atomic species. The PBE version of the generalized 
gradient approximation (GGA)^^ is used for the elec- 
tron exchange and correlation, and optimized TrouUier- 
Martins pseudopotentials^ are used for the atomic cores. 
The atomic structure of the molecule, lead surfaces, and 
molecule-lead separation are fully optimized, with resid- 
ual forces less than 0.02 eV/A. 

For the transport calculationSSi we divide an infi- 
nite LML system into three parts: left lead (L), right 
lead (R), and device region (C) which contains the 
molecule and large parts of the left and right leads, 
as shown in Fig. 1, so that the molecule-lead interac- 
tions can be fully accommodated. For a steady state 
situation in which the region C is under a bias Vb 
(zero or finite), its density matrix (Dc) and Hamilto- 
nian (He) can be determined selfconsistently by the 
DFT+NEGF techniques.iLi2*i2i2a-2ii2^ The Kohn-Sham 
wave-functions are used to construct a single-particle 
Green function from which the transmission coefficient 
at any energy, T{E, Vb), is calculated. The conductance, 
G, then follows from a Landauer-type relation. 



III. RESULTS AND DISCUSSION 
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FIG. 2: Comparison among the transmission functions of the 
S-, Se-, and Te-anchored systems: (a) (001) with adsorption 
at the hollow site, (b) (111) hollow site, (c) (111) bridge, and 
(d) (111) top site. The conductance is largest in the (001) 
case, and decreases as the atomic number of the anchoring 
group increases. 



A. Structure 

We find that if the anchoring atom is adsorbed at the 
hollow site, the effect of contact atomic relaxation is mi- 
nor and independent of the lead orientation [see Fig. 1 
(a) and (b)]. However, if the anchoring atom is adsorbed 



at the b or t site [for a Au(lll) lead], the relaxation effect 
becomes significant [Fig. 1 (c) and (d)] . This behavior is 
understandable because the h site is the bulk atomic po- 
sition, so that the adsorption of an anchoring atom at it 
will not noticeably change the directional binding charac- 
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ter of the surface. For the b or t site, in contrast, the di- 
rectional binding character can be significantly modified 
by an adsorbed atom, leading to the significant contact 
atomic relaxation shown. 

The relative stability of the diff'erent adsorption sites is 
related to the binding energy between the molecule and 
the lead, given by the difl^erence in total energy between 
the whole system and the two subsystems using the same 
supercell and k-sampling. Table I summarizes the results 
for both (001) and (111) leads connected by the three an- 
choring atoms, both relaxed and unrelaxed. We wish to 
point out three features: First, the binding energy on 
(111) is larger than that for (DDI). This is a straight for- 
ward result of the different directional binding character 
of the two surfaces with the group- VI anchoring atoms: 
the coordination on the (111) surface is smaller than on 
(001) and so is more favorable to the group- VI anchor- 
ing atom. Second, for the Au(OOl) lead, the most stable 
adsorption site, regardless of contact relaxation, is the h 
site [shown in Fig. 1 (a)] followed by b and then t. The 
energy gains for b and t are quite large: for the relaxed 
S-anchored system, for instance, the energy increase for 
b and t sites are 0.77 and 1.80 eV/contact, respectively. 
Third, for the Au(lll) lead the situation is not so simple. 
If the contact is not relaxed, the most stable adsorption 
site is h followed by b and t; however, after relaxation, the 
most stable site becomes b followed by h and t. The en- 
ergy differences are smaller here compared to the (001) 
case: for the unrelaxed S-anchored system, the energy 
difference between h and b is 0.41 eV/contact while after 
relaxation it is only —0.21 eV/contact. 

Our results are largely in agreement with previous re- 
sults concerning chemisorption on Au(lll). Short alkane 
thiols on unrelaxed Au(lll) prefer the h sitei^ In a study 
of the chemisorption of 1,4-benzene-dithiol perpendicu- 
larly on the surface of an unrelaxed Au25 cluster;^ the 
most stable adsorption site was h followed by b then t, as 
in our results. Finally, a systematic ab initio calculation 
of S-CeHs chemisorbed on a AU29 clusterjSi in which the 
whole structure was fully relaxed in different geometrical 
configurations, it was found that several b sites have the 
lowest energy, consistent with our result for the relaxed 
S-anchored system. Despite this agreement with previous 
results, we caution that the structural optimization and 
energetics of the different structures in this paper may 
include a significant contribution from the small width 
of the Au leads; as a result, the present results may dif- 
fer quantitatively from those for chemisorption on either 
an infinite surface or a small Au cluster. 

Considering that in break-junction experiments differ- 
ent contact atomic structures will occur and an adsorbed 
molecule may not be at its global equilibrium position, we 
investigate the three adsorption sites on Au(lll)-leads. 
For Au(GGl)-leads, however, we only investigate the most 
stable adsorption site (h) because it has a substantially 
lower energy than the next most stable one. 



TABLE II: Equilibrium conductance (G, in units of 2e^//i) 
and molecule-lead electron transfer (AQ, in units of electron 
charge). A positive AQ means that electrons are transferred 
from lead to molecule. 'h-Au' inidcates that the molecule is 
connected to the leads via an additional Au atom adsorbed 
at the hollow site. The sensitivity to both atomic structure 
and type of anchoring atom shows that treating the contact 
realistically is critical for reliable calculations in molecular 
electronics. 



anchoring 


lead 


adsorption 






atom 


orientation 


site 


AQ 


G 


s 




h 


-0.048 


0.053 




fill) 


h 


+0.053 


0.008 






b 


+0.117 


0.002 






t 


+0.237 


0.013 






h-Au 


+0.228 


0.490 




jellium 




+o.r 


0.036" 


Se 


(001) 


h 


-0.206 


0.031 




(111) 


h 


-0.010 


0.005 






b 


+0.036 


0.002 






t 


+0.170 


0.009 






h-Au 


+0.235 


0.550 




jellium 






0.12° 


Te 


(001) 


h 


-0.294 


0.014 




(111) 


h 


-0.044 


0.002 






b 


-0.002 


0.001 






t 


+0.167 


0.004 






h-Au 


+0.240 


0.430 




jellium 




-0.5° 


0.88° 



Result from Refii^ 



B. Effects of structure on transmission 

In Fig. 2 we show for the systems studied the trans- 
mission function T{E) under zero bias. The values of 
the equilibrium conductance are summarized in Table II, 
together with the molecule-lead electron transfer deter- 
mined by a MuUiken population analysis. 

In Table II, first note that the Au(OOl) lead gives a 
much larger conductance (by about 6 times) than the 
(111) lead for all three anchoring atoms - S, Se, and Te. 
This can be understood by analyzing the molecule-lead 
coupling and electron transfer for the two lead orienta- 
tions: In the (001) case, the anchoring atom has four 
nearest neighour Au atoms while there are only three 
for (111). This difference in the contact atomic configu- 
ration certainly affects the molecule-lead coupling, which 
may also lead to a difference in the molecule-lead electron 
transfer. However, for the S-anchored systems, both sur- 
faces result in very small electron transfer (Table II) , but 
their transmission functions are essentially different [Fig. 
2 (a) and (c)]. In particular, there is a small peak around 
the Fermi energy for (001) which is absent for (111), in- 
dicating that the difference in equilibrium conductance is 
mainly due to the different molecule-lead couplings. An 
interesting aspect here is the difference between the bind- 
ing (Table II) and the transport with the Au(OOl) and 
Au(lll) leads: contacts with good (poor) binding have. 
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however, relatively poor (good) transparency. This is an 
illustration that for completely different contact struc- 
tures, binding and transport character are not necessarily 
related. 

Another difference between the two lead orientations is 
that the overall structure of T{E) is totally different (Fig. 
2): the T{E) for Au(OOl) are quite smooth while those 
for Au(lll) have many sharp features. Possible reasons 
for the latter include a Au(lll) lead that is too thin or 
its much lower symmetry. Note that other calculations^^ 
also find sharp structure when using a Au(lll) lead. To 
dramatize the difference in transport between the two 
lead orientations, we calculate T{E) for infinite pure 
(001) and (111) wires (Fig. 3). Because of the infinite 
periodic structures of the two wires, both T{E) are step 
functions. As can be seen, the average transmission of the 
(001) wire is larger than that of the (111) wire, and the 
(111) transmission coefficient fiuctuates strongly. These 
strong fluctuations are related to the sharp features in 
the transmission through the molecule. 

Table II shows a clear trend in the equilibrium conduc- 
tance for the different adsorption sites: G{t) > G(h) > 
G(b) for all the S-, Se-, and Te-anchored (111) systems. 
Comparing to the binding energies for the different sites 
(Table I), we see again that there is no direct relation be- 
tween transport and binding properties when comparing 
different structures. Note that the maximum difference 
in G for the three sites can be up to about 6 times. 

In contrast, the trend in molecule-lead electron transfer 
is AQ{t) > AQ{h) > AQ{h). Here we can also find cases 
where the electron transfer for different adsorption sites 
is very close but the conductance is quite different [in Ta- 
ble II compare Se_(lll)_h to Se_(lll)_b, or Te_(lll)_h to 
Te_(lll)_b]. These indicate that the trend in the equilib- 
rium conductance is mainly due to differing molecule-lead 
coupling for the three adsorption sites rather than simply 
differing charge transfer. 

This conclusion can also be reached from Fig. 2. In 
Fig. 2 (b) and (c), because of the similar electron trans- 
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E-Ef (eV) 

FIG. 3: Transmission functions of infinitely long wires: 
Au(001)-2V2x2V2 (solid line) and Au(lll)-2x2 (thin line). 
Note the much larger average transmission for (001) and the 
large fluctuation for (111). 



TABLE III: Current (in ^A) for several systems under 1 V 
and 3 V bias. The notation is as in Table II. 



anchoring 


lead 


adsorp. 




current 


atom 


orient. 


site 


1 V 


3 V 


S 


(001) 


h 


6.71 


44.82 




(111) 


b 


1.57 


20.77 






h-Au 


2.29 


7.36 


Se 


(001) 


h 


4.74 


32.57 




(111) 


b 


1.29 


14.81 






h-Au 


2.37 


6.40 


Te 


(001) 


h 


2.37 


17.43 




(111) 


b 


0.94 


9.13 






h-Au 


2.00 


5.13 



fer (Table II), T{E) for the h- and b-adsorbed systems 
are similar around the Fermi energy. In addition, the 
relative positions of the Fermi energy in the gap (which 
is comparable to the ~ 0.3 eV HOMO-LUMO gap of the 
molecule including the two anchoring atoms) are similar. 
However, there is a small shoulder at the Fermi energy 
for the h-adsorbed cases which is absent in the case of 
b-adsorption; this is the signature of the difference in the 
molecule-lead coupling. 

Results for finite bias provide further support for these 
conclusions. In Table III we list the current under 1 and 
3 V biases for the (001)_h and (lll)_b systems with each 
of the anchoring atoms. The first thing we should notice 
is that the amplitude of the current for the S-anchored 
system is consistent with both a previous selfconsistent 
DFT-t-NEGF calculation using a cluster geometry2SiS& 
and a self-consistent jellium model calculationsi^ All 
of these theoretical results are much larger than ex- 
perimental rcsults'^^ (as well as the result from a non- 
selfconsistent DFT+NEGF calculation) /^S, The reason for 
this discrepancy between theory and experiment is still 
an open problem; it may be related to the structural 
difference between the theoretical model and the real ex- 
perimental conditions. 

Under finite bias the current for (001) is much larger 
than that for (111), although the relative difference tends 
to decrease with increasing bias. To show directly differ- 
ences between the two lead orientations under a finite 
bias, the voltage drops for 1 V bias between the two lead 
surfaces are shown in Fig. 4 for the (001)_h and (lll)_b 
cases. Despite the difference in conductance between the 
two lead orientations, the voltage drops are strikingly 
similar. Note features in the middle corresponding to 
the benzene ring, features between the anchoring atom 
and the surface, and features between the anchoring atom 
and the ring. Looking carefully in Fig. 4 [(a) vs. (d), 
(b) vs. (e), and (d) vs. (f)], we see that the voltage 
drop between the left surface and the anchoring atom is 
slightly sharper in the (111) systems than in the (001) 
cases, indicating that the surface-anchoring-atom cou- 
pling is slightly weaker at the (111) contact than at the 
(001) contact. 
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FIG. 4: Voltage drops for 1 V bias (between the two lead surfaces as shown in Fig. 1). The contour surfaces are plotted in the 
plane crossing the benzene ring. The first column is for (001)_h cases: (a) S_(001)_h, (b) Se_(001)_h, and (c) Te_(001)_h. The 
second column is for (lll)_b cases: (d) S_(lll)_b, (e) Se_(lll)_b, and (f) Te_(lll)_b. Note the similar voltage drop behavior 
for the different anchoring atoms and lead orientations. 



C. Chemical trends 

The chemical trends of different anchoring groups is a 
particularly interesting and important aspect of molec- 
ular electronics They demonstrate, for instance, the 
critical role of realistic contact atomic structure in elec- 
tron transport through molecules. In Fig. 2 and more 
clearly in Table II, a chemical trend is evident in both 
the molecule-lead electron transfer and the equilibrium 



conductance for the three anchoring atoms we study, re- 
gardless of the lead orientation or adsorption site: As the 
atomic number of the anchoring atom increases, S — > Se 
Te, both the electron transfer from the leads to the 
molecule and the conductance decrease. 

Our result of decreasing G with increasing atomic num- 
ber is directly opposite to the conclusion reached by a 
previous calculation^^ using the jellium model for the Au 
leads (Table II). Besides the qualitative difference in the 
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TABLE IV: Equilibrium conductance (in units of 2e^/ft) for 
the unrelaxed systems with larger Au(OOl) leads: 3\/2 x 3\/2, 
4v^ X 4a/2, and 4v^ x 4\/2 periodic surface. Note the same 
chemical trend for the three anchoring atomsas in the case of 
the smaller 2^2x2^ leads (Table II). 



anchoring 




lead size 




atom 


3\/2 X 3%/2 


4^2 X 4^2 4^2 


X 4\/2 surface 


S 


0.110 


0.089 


0.108 


Se 


0.076 


0.064 


0.087 


Te 


0.044 


0.037 


0.053 



chemical trend of conductance, the present calculation 
shows that the maximum difference in equilibrium con- 
ductance is less than 5 times in contrast with 25 times 
found in the previously. The chemical trend in the elec- 
tron transfer is consistent with the electronegativities of 
the three species [2.58, 2.55, and 2.1, respectively (Paul- 
ing scale)]; it is also consistent with the trend in binding 
energy (Table I) , indicating that for similar contact struc- 
tures the transport and binding properties are related. 

In order to check whether the clear chemical trend in 
G is affected by the width of the leads, we carry out 
calculations for (001)_h systems using two wider leads: 
3V2 X 3V2 and 4\/2 x 4\/2. Furthermore, we also carry 
out test calculations using a (001)-(4V2 X 4^2) periodic 
surface for the leads, which will be a good approxima- 
tion to an infinitely large surface because of the large 
separation between the molecule and its images (larger 
than I2A). To avoid prohibitive computational effort, we 
adopted a single zeta plus polarization (SZP) basis set 
here. Table IV shows that the chemical trend remains 
with the wider Au leads. 

Results for finite bias provide further support for our 
conclusions. Under 1 or 3 V bias, the current decreases 
with increasing atomic number of the anchoring atom 
for both (001) and (111) leads (Table III). T{E), shown 
in Fig. 5 for the (001)_h cases, decreases at most ener- 
gies as atomic number increases. From the voltage drop 
profiles in Fig. 4, for both (001) and (111) the drops 
around the anchoring atoms become slightly sharper as 
the atomic number of the anchoring atom increases. This 
is direct evidence that the surface to anchoring-atom cou- 
pling under finite bias becomes weaker as atomic number 
increases. 

In break junction experiments, the break surfaces are 
certainly not as flat as assumed above but rather very 
likely have atomic fluctuation (atomic scale roughness). 
To see whether or not the clear chemical trend of the dif- 
ferent anchoring atoms is affected by the change in con- 
tact atomic configuration, we consider an additional Au 
atom adsorbed at the h site of the two Au(lll) lead sur- 
faces. There has been previous work on this system us- 
ing DFT-I-NEGF calculation32S*2Si2i using an unrelaxed 
S_(lll)_h configuration, as well as a calculation model- 
ing the leads with jellium^ First, using an approximate 
non-self-consistent approach in which a constant imagi- 
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FIG. 5: Comparison among the transmission functions for 
S-, Se-, and Te-anchoring in the case of (001) adsorption at 
the hollow site: (a) 1 V bias, (b) 3 V bias. Note that the 
transmission coefficient decreases as the atomic number of 
the anchoring group increases, the same trend as for zero bias 
(see Fig. 2). 



nary self-energy is adopted for the Au lead, one study 
found that the current under 4 V bias depends strongly 
on the adsorption site on the Au lead''^; when the con- 
nection was made by an apex Au atom, the current was 
significantly reduced. Second, in a systematic calculation 
for the unrelaxed S_(lll)_h system adopting a cluster 
method in which only 6 Au atoms are included to form 
the device regiou j^^i^^ the equilibrium conductance was 
significantly increased by the addition of the Au atoms 
at each contact, but the current under large bias was 
decreased. The result for equilibrium conductance was 
opposite to that obtained with jellium leadsji^ showing 
that one must use caution when considering the latter. 

In view of these previous results, our goal here is two 
fold: On the one hand, we wish to study this system with 
the more recently developed presumably more accurate 
methods; in particular, we will include 36 Au atoms per 
(111) lead and 50 per (001) lead in contrast to the 6 used 
previously. Summarizing the results below, we find that 
the more accurate methods give qualitatively the same 
results as those used previously. 

On the other hand, and more importantly, we wish to 
investigate the chemical trend in the effect of the addi- 
tional Au atom. Our results for electron transfer and 
equilibrium conductance are listed in Table II, with the 
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FIG. 6: Voltage drops for 1 V bias with an additional Au 

atom at each contact in the h site on (111): (a) S (b) S.and 
(c) Te. The additional Au atom is included as a part of the 
molecule not as a part of the surface. The contour surfaces 
arc plotted in the plane crossing the benzene ring. Note that 
the voltage drop behavior is very similar for the three different 
anchoring atoms but very different from that in the systems 
without the additional Au atoms (see Fig. 4). 



adsorption site denoted 'h-Au'. The additional Au atom 

increases significantly the electron transfer from the leads 
to the molecule, while the difference among the three an- 



choring atoms is strongly reduced. Because of the large 
electron transfer, the LUMO resonance becomes close to 
the Fermi energy; consequently, the equilibrium conduc- 
tance also increases significantly. Note that introduction 
of the additional Au atoms at the contacts changes the 
chemical trend: Se becomes the best anchoring atom in 
terms of contact transparency. 

Under a finite bias, the relatively large amount of 
charge transferred to the molecule causes a large den- 
sity of states at energies between the left and right Fermi 
energies. This in turn causes very different voltage drop 
behavior in the presence of additional Au atoms: fluctu- 
ations in the voltage drop profile become large because 
of a large bias-induced polarization (compare Fig. 6 with 
Fig. 4). If the bias is relatively small (1 V in Table III), 
Se remains the best anchoring atom for contact trans- 
parency; however, if the bias is large (3 V), then the 
chemical trend reverts to S being best. Although the 
introduction of the additional Au atoms increases signif- 
icantly the molecular conductance for zero or small bias, 
it decreases the current under large bias (Table III). 



IV. SUMMARY 

By using a state-of-the-art ab initio method for elec- 
tronic structure and electron transport, we have carried 
out a systematic calculation for the molecular conduc- 
tance of benzene sandwiched between two Au electrodes. 
Our calculation is the first to fully include the effects 
of relaxing the contact atomic structure. The main re- 
sults are: (1) Detailed contact structure strongly affects 
molecule-lead coupling, electron transfer, and molecu- 
lar conductance. (2) There is no general relation be- 
tween the binding strength of an anchoring atom and the 
transparency of the contact; however, in the special case 
of comparing contacts with very similar atomic struc- 
ture, stronger binding does imply increased transparency. 
(3) For ideally fiat break surfaces, the equilibrium con- 
ductance decreases with increasing atomic number of the 
anchoring group regardless of the adsorption site, lead 
orientation, or bias. (4) This chemical trend is, however, 
affected by the local contact atomic configuration: An 
additional Au atom at the contact with the Au(lll) lead 
changes the best anchoring atom for small bias (from S 
to Se) although not for large bias. These results demon- 
strate the critical role of realistic contact atomic structure 
in electron transport through molecules. 
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